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We investigate degree of coherence of pion sources produced in relativistic heavy ion collisions using multi- 
particle interferometry. In order to obtain "true" chaoticity, X tme from two-pion correlation functions measured 
in experiments, we make a correction for long-lived resonance decay contributions. Using this X tme and the 
weight factor which are obtained from parameter fitted to two- and three-pion correlation function, we cal- 
culate a chaotic fraction e and the number of coherent sources a for different colliding energies. The result 
gives constraints on the source and shows an increase of the minimum number of a with multiplicity, although 
multiplicity independent chaoticity is not excluded. 
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I. INTRODUCTION 

Relativistic heavy ion collisions provide us an unique op- 
portunity to explore nature of hot and dense nuclear matter on 
the earth. In the highest energy collisions at the BNL-RHIC, 
it is expected that the matter created soon after the collision 
of two heavy nuclei can be the strongly interacting quark- 
gluon plasma, which gradually cools down and then becomes 
hadronic matter via phase transitions. To understand the na- 
ture of the QCD matter, it is important to know what kind of 
information experimental observables contain. Pion interfer- 
ometry has been one of the most important observables be- 
cause it can give us information on sizes of the source which 
pions come from, through the HBT effect. The HBT effect is 
a quantum-mechanical effect due to symmetrization of two- 
boson wave function and occurs if the source is not completely 
coherent. The strength of the two-boson momentum intensity 
correlation takes its maximum value in the case of the per- 
fectly chaotic source. Because the chaoticity can be related to 
the degree of thermal equilibrium and to composition of the 
source, analyses of the chaoticity can provide information on 
a state of hadronic matter which may reflect how hadrons are 
produced. 

Experimentally measured two-pion chaoticity 
A, = C2(p,p) — 1 does not usually reflect real coherence 
of the source because of long-lived resonance decay contri- 
butions. Three-pion correlation function is more useful for 
this purpose. The three-pion correlator measured in some 
experiments is 

r 3 (<23) = 

[C 3 (6 3 ) - 1] - [C 2 (Qn) - 1] - [C 2 (g 23 ) - 1] - jPzgfei) - 1 
VlCiiQn) - l][C 2 (e 2 3) - l][C 2 (63i) - 1] 

(1) 

where Q y = ^-(pi-pj) 2 , Q 3 = y 1 Q 2 U + Q\ x , and C„ 
denotes the n-particle correlation function. The weight factor 
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CO = rj (0)/2 characterizes the degree to how pions sources are 
chaotic. For a completely chaotic source, co = 1. If the source 
chaoticity is really a physical quantity, both two- and three- 
pion correlation function should give quantitatively consistent 
value of the chaoticity. However, experimental results do not 
seem so, because of the apparent reduction of the X due to 
long-lived resonances. But it has been shown that we can im- 
pose stronger constraints by using both of two- and three-pion 
correlation data 0]] . 

In this work, we find that such decay contributions can be 
eliminated by making use of a statistical model [2] and applied 
it to various experimental data measured in SPS and RHIC ex- 
periments yfl . Then, we obtain the weight factor co from two- 
and three-pion correlation function data in these experiments. 
Using the resonance-corrected A, which we call A tme and CO, 
we calculate a chaotic fraction e and the number of coherent 
sources based on a partially multicoherent source model |4]. 
From this result, we discuss how the structure of the sources 
changes from low energy collisions at SPS and higher ones at 
RHICJI. 



n. CORRECTION OF X FACTOR FOR LONG-LIVED 
RESONANCE DECAY 

In the presence of long-lived resonance decay contributions 
to two-pion correlation function, the two-pion chaoticity for a 
chaotic source is reduced as 

/ N 1 

V^=l- — (2) 

where N K is the total number of pions and 7V^ is the number 
of pions coming from long-lived resonance decay [6]. Hence, 
if one can estimate the ratio of pions from the decay to the 
total number, one can obtain apparent reduction factor caused 
by the decay contribution. Indeed, it is possible with the help 
of the statistical model which has been known to give a nice 
description of observed particle number ratio. Under the as- 
sumption of the constant chemical freeze-out temperature and 
chemical potential and boost-invariance along collision axis, 
the ratio of particle number can be written as the ratio of num- 
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ber density, i.e.,Ni/Nj = ni/rij, where is given by 

n i = ^J"dpp 2 f(E,T,n) (3) 

with f(E,T,fij being the equilibrium distribution functions 
|0]. Using these formulae, we calculate the ^ eff for S+Pb 
and Pb+Pb collisions at the SPS and Au+Au collisions at 
the RHIC and obtain A, tme through a relation A. tlue = X ex P/l eff 
where X exp is momentum-averaged experimental data Q,0]. 
The temperature and baryonic chemical potential are deter- 
mined by the y} fitting to experimental data of particle ratio. 
See Ref. 10] for details. Here we summarize the result in Ta- 
ble m 

TABLE I: Summary of l tme 



System 


(7> B )[MeV] 




^ true 


S+Pb (NA44) 


(173,196) 


0.59(6) 


0.94(6) 


Pb+Pb (NA44) 


(161,223) 


0.55(3) 


0.98(3) 


Pb+Pb (WA98) 


(161,223) 


0.58(4) 


1.03(4) 


Au+Au (STAR) 


(158,36) 


0.57(6) 


0.93(8) 



III. EXTRACTION OF THE WEIGHT FACTOR 

In order to obtain the weight factor co, we have to extrapo- 
late experimentally measured rj(Qj,)[Eq. (fl~|i] to Qi = 0. Us- 
ing a simple source function in which instantaneous emission 
and spherically symmetric source are assumed, we construct 
the two- and the three-pion correlation functions with the for- 
mulae 



Ci{Pl,P2,P"i) 




(4) 



(5) 



where fij = 1/y cosh(7?|q l; |) is the source function and R, 

^i nv and V are parameters which should be determined by a 
simultaneous % 2 fit to the two- and the three-pion correlation 
functions ill Oil . From a set of these parameters, we can calcu- 
late CO using Eqs. dUi,© and (Q}. The results for co are shown 
in Table HH 



we investigate the coherence of the sources using the partially 
multicoherent model [4]. In this model, there are two charac- 
teristic parameters which are related to the A. true and co as 



^true 



CO : 



a 



a+(l-e) 2 ' 

2a 2 + 2a(l -e) 2 + 3(l -e) 3 (l -2e) 
2[a 2 + 3a(l -e) 2 + (l -e) 3 ] 



(6) 



a- 



(1-e) 2 



a 



(7) 



By solving the above equations with respect to e and a (this 
can be analytically), we can obtain allowed regions for e and 
a corresponding to the available range of A, tme and CO. 




n 0.2 0.4- 06 0.8 



FIG. 1 : Allowed region for e and a. From top to bottom, each figure 
(a)-(d) shows S+Pb, Pb+Pb by NA44, Pb+Pb by WA98 and Au+Au 
case, respectively. 



TABLE II: Results for CO 



S+Pb(NA44) Pb+Pb(NA44) Pb+Pb(WA98) Au+Au(STAR) 



co 0.40±0.44 1.15±0.67 



0.78±0.44 0.872±0.097 



IV. ANALYSIS WITH THE PARTIALLY 
MULTICOHERENT MODEL 

From considerations in previous sections, we have obtained 
the two quantities, A. tme and CO as experimental results. Next, 



The result is shown in Fig. Q] In each of figures, the light- 
est shade area labelled "A" and the second one labelled "B" 
denote the allowed regions coming from X tme and CO, respec- 
tively. The darkest areas labelled "C" are overlap of area "A" 
and "B", then correspond to the allowed parameter region for 
e and a. The best fit points are indicated by the filled box. 

From Fig. Q] it seems to be difficult to find systematic 
change of the allowed regions. This result mainly comes 
from the fact that A. tlue 's are close to unity in the Pb+Pb data. 
However, it has been suggested that Coulomb correction to 
the two-pion correlation functions is over-corrected one; the 
values of A, exp can decrease if we take account of the partial 
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Coulomb correction [8, 9]. Hence, if appropriate corrections 
were made, obtained A. tlue becomes smaller. In order to obtain 
a rough sketch of a tendency, we also draw the allowed regions 
for A,"™ which is 0.8 times smaller than A, tme . For S+Pb data, 
we multiply the original A, tme by factor 0.7 because Gamow 
(point-like source) correction is made for this data. 



l lrur =Q.99B>[).7' 
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FIG. 2: Allowed region for £ and a, in the case of using From 
top to bottom, each figure (a)-(d) shows S+Pb, Pb+Pb by NA44, 
Pb+Pb by WA98 and Au+Au case, respectively. 

From Fig. |2 we can see that all allowed regions (the darkest 
shaded areas) are narrow but there seem to exist systematics. 
The best fit point seems to move upper left (small chaotic frac- 
tion and large number of coherent sources) side, except for the 



NA44 Pb+Pb result. Most important feature of this result is 
that the upper limit of e and lower limit of a are determined 
by the lower limit of CO. We plot the maximum and the min- 
imum values of a in Fig. [3] The clear increase of minimum 
number of the coherent sources can be seen as a function of 
multiplicity while maximum number of those shows no such 
clear tendency. 

In summary, we have given an analysis of the degree how 
chaotic the pion sources are in relativistic heavy ion collisions 
at the SPS and the RHIC. The analysis can be done by using 
both two-pion correlations and three-pion correlations. We 
find that the correction for long-lived resonance decay contri- 
butions to the two-pion correlation function can be subtracted 
with the help of the statistical model. From a point of view in 
which multicoherent sources and a background chaotic source 
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FIG. 3: Maximum and minimum values of a as a function of multi- 
plicity. For Pb+Pb data, both plots are slightly shifted along horizon- 
tal axis for a clear comparison of the results. 



are produced, we show that the model gives constraints on the 
structure of the source. Although the maximum number of 
the coherent sources does not show a clear multiplicity de- 
pendence, the minimum number of coherent source increases 
as the multiplicity increases. 
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